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Background of the Invention 
[0001] 1. Field of the Invention. 

[0002] This invention generally relates to acoustic waveguides and in particular to a 
method and system for modeling the design of an acoustic waveguide based upon 
predicted performance standards and performance metrics for a waveguide having certain 
physical characteristics and dimensions. 

[0003] 2. Related Art. 

[0004] Often times, loudspeakers consist of a transducer or driver unit coupled to a 
waveguide. A waveguide can also be commonly referred to as a horn or acoustic 
waveguide. A waveguide functions to provide gain for the transducer, i.e., increases the 



acoustic sensitivity of the loudspeaker in a region of frequencies. A waveguide can also 
assist in the control of dispersion on and off-axis as well as assist with directivity mating 
with other transducers and can simplify loudspeaker system integration. 
[0005] Typical waveguides include a "throat" or entrance at one end and a "mouth" 
or exit at the opposing end. The throat end of the waveguide is typically coupled to the 
transducer or driver and receives the initial input of sound from the driver. The 
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waveguide then usually increases in cross-sectional area or flares out as it approaches the 
mouth. The sound is then dispersed through the mouth, which is the exit of the 
waveguide. Thus, the throat end of the waveguide is typically narrower in cross-section 
in both the horizontal and vertical directions and generally defines a bounded region that 
directs the sound from the throat to the mouth of the waveguide. This interior bounded 
region may be referred to as the waveguide profile. The sound produced as planar 
surfaces parallel to the throat, are referred to as wave fronts. 

[0006] In operation, the surfaces of the waveguide in a loudspeaker typically produce 
a coverage pattern of a specified total coverage angle that may differ horizontally and 
vertically. T he c overage angle i s a total angle in any p lane o f observations, although 
horizontal and vertical orthogonal planes are typically used. The coverage angle is 
evaluated as a function of frequency and is defined to be the angle at which the intensity 
of sound (Sound Pressure Level - SPL) is half of the SPL on the reference axis, which is 
the axis direction usually normal to the throat of the driver. 

[0007] Acoustic energy radiates into the throat from the transducer at high pressure, 
with a wave front that is nominally flat and free of curvature. As the wave front expands 
outward to toward the mouth of the waveguide, the axial area increases in a uniform and 
monotonically increasing fashion. Analogous to electrical transformers in the electrical 
domain, waveguides can be considered as acoustical transformers in the acoustical 
domain. In the acoustical domain, waveguides contain impedance along the profile with 
resistive and reactive components. However, sound pressure level is produced primarily 
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by the acoustical resistance of the waveguide. That is, acoustical reactance does not 
contribute to the sound pressure level In the work presented, the rate of increasing area is 
controlled by an area expansion function designed to provide minimal acoustic reactance 
(or m aximum a coustic r adiation r esistance a 1 1 he t hroat). This a pproach i ncreases t he 
sensitivity and ultimately, the efficiency of the transducer and waveguide assembly. 
[0008] The determined area expansion rate is intended to create a uniform dispersion 
pattern on and off-axis by manipulating the acoustical impedance as a n function of 
frequency to theoretically lower frequency range of operation. The coupling of the 
waveguide acoustic impedance source to the acoustic impedance of the surrounding 
environment; provides an action analogous to an electrical transformer. The winding 
ratio is equivalent to the ratio of the radiation resistance seen by the driver and the 
radiation resistance of the surrounding environment. In this analogy, the change in 
pressure from the throat to the mouth of the waveguide is equivalent to the change in 
voltage across an electrical transformer. 

[0009] The shape of an acoustic waveguide affects the frequency response, polar 
pattern and the level of harmonic distortion of sound waves as they propagate away from 
the acoustic waveguide. As loudspeakers produce sound waves, waveguides are used to 
control the characteristics of the acoustic wave propagation. As previously stated, the 
increase in area of the waveguide from throat to mouth is typically controlled by an area 
expansion function designed to provide appropriate acoustic impedance. Many different 
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theories on waveguide design have been developed in the past to help determine the 
optimal expansion functions for waveguide designs. 

[0010] One c ommon d esign approach, d eveloped b y K eele, i nvolves a t wo-section 
waveguide or horn design. In this design approach, an exponential design is used on the 
section near the throat, while the outer section utilizes a conical design approach. 
Similarly, Geddes developed an alternative design approach that is a well known in the 
industry. This approach uses exponential algebraic equations and functions developed by 
Geddes to determine the optimal contour of a waveguide once required v alues for the 
throat radius and coverage angle have been determined. 

[0011] Current design approaches, such as those taught by Keele and Geddes, first 
determine the desired performance standards of the waveguide and then design the 
waveguide using established exponential functions or algebraic equations that are 
designed to model a waveguide to achieve the desired standards. No design method 
currently exists, however, that uses the performance standards of a waveguide of known 
contours and dimensions as a design metric. Additionally, no design method currently 
exists that captures the change in acoustic impedance, in particular the change in acoustic 
reactance, along the profile of the waveguide as part of the design standard. A need 
therefore exists for a waveguide design method such that one can predict the performance 
standards of waveguides having various contours and dimensions without the necessity of 
building a prototype. Under this proposed approach, design iterations can be made 



4 



U.S. Express Mail No. EU 1 24349650 US PATENT 
Filing Date: October 29, 2003 Attorney Docket No. HI03025USU 

(P04009USU) 

before the prototype stage of the waveguide since the performance standards may be 



predicted in advance of the design. 
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Summary 



[0012] This invention provides a method of designing waveguides capable of 
sustaining a generally constant change in impedance and pressure gradient along the 
transition of the waveguide from throat to mouth by using design metrics known to 
correlate with the physical dimensions, contours, and acoustical measurements of 
waveguides. The design methodology captures the change in acoustical impedance 
within the area expansion function and explicitly determines the waveguide profile 
required by providing a predicted frequency response, without the use of a discrete 
prototype. 

[0013] With an established set of design metrics, waveguide profiles can be design by 
dividing the waveguide profile into two or more different exponential profiles having two 
or more different slopes. The slopes are then altered by applying functions derived from 
the set of design metrics. Once altered, the resulting waveguide profiles from the 
different slopes are then concatenated together and smoothed to produce a design key for 
prototyping a waveguide that can achieve the desired design performance specifications; 
for which the design metric is based. 

[0014] In one embodiment, the design metric is the change in acoustic reactance 
along the profile of the waveguide. The waveguide is divided into ten sections. Initial 
values are then assigned for the radius or diameter of the throat of the waveguide as well 
as values for the initial slope of the waveguide along the major and minor (or x and y) 
axis, polynomial smoothing order for the ten concatenated profiles, and the desired depth 
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of the waveguide. The values for the slopes of each section are then altered based upon 
functions derived from the design metrics. In this example implementation, each slope is 
adjusted to minimize the change in acoustic reactance along the waveguide profile, which 
is the desired performance standard. Once the slopes of each section are adjusted to 
achieve minimal change in acoustic reactance, the sections are concatenated together and 
the curve is smoothed using a polynomial function order curve fit to create a continuous 
waveguide profile. The profile correlates with the design measurements, which allows 
for the prediction of the performance standards or dispersion characteristics of the 
waveguide. Design iterations may then be made to adjust for desired performance 
measurements without the necessity of building a prototype. 

[0015] Furthermore, since the uniform acoustical reactance along the waveguide 
profile provides stable and predictable dispersion on-axis and off axis, the invention may 
be used to design waveguides having elliptical cross-sectional areas that produce circular 
dispersion patterns (i.e. an elliptical waveguide that produces the same horizontal and 
vertical dispersion patterns from 1kHz to 10kHz). Conversely, the design allows for the 
design of waveguides having circular cross-sectional areas yet provide elliptical 
dispersion patterns (i.e. a circular waveguide that produces different horizontal and 
vertical dispersion patterns from 1kHz to 10kHz). 

[0016] Other systems, methods, features and advantages of the invention will be or 
will become apparent to one with skill in the art upon examination of the following 
figures and detailed description. For example, this design method could be used to 
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design transducers diaphragms found in tweeters, mid-ranges, mid-bass, woofers, and 
subwoofers commonly used in loudspeaker systems. Similarly, this work could be used to 
design waveguides that are found in radar and communication applications using 
analogous partitions, concatenations, and design metrics. It is intended that all such 
additional systems, methods, features and advantages be included within this description, 
be within the scope of the invention, and be protected by the accompanying claims. 
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Brief Description of the Figures 



[0017] The invention can be better understood with reference to the following 
figures. The components in the figures are not necessarily to scale, emphasis instead 
being placed upon illustrating the principles of the invention. Moreover, in the figures, 
like reference numerals designate corresponding parts throughout the different views. 
[0018] FIG. 1 is a front view of a loudspeaker utilizing an acoustic waveguide 
designed in accordance with the design method of the invention. 

[0019] FIG. 2 is a cross-sectional view of the acoustic waveguide and dome 
diaphragm of the loudspeaker of FIG. 1 taken along line A — A. 

[0020] FIG. 3 is a flow diagram of an example implementation of the waveguide 
design methodology of the invention. 

[0021] FIG. 4 illustrates the transition in theoretical component of the acoustic 

impedance along the transition of a waveguide from throat to mouth. 

[0022] FIG. 5 illustrates the x-axis, y-axis, z-axis and the radius r for an example 

waveguide. 

[0023] FIG. 6 illustrates a waveguide divided into ten (10) sections. 

[0024] FIG. 7 illustrates a depth verses height profile of an example waveguide 

designed in accordance with the invention. 

[0025] FIG. 8 illustrates the predicted change in acoustic reactance verses frequency 
along the x and y axis for the waveguide profile of FIG. 7. 
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[0026] FIG. 9 illustrates a depth verses height profile of another example waveguide 
designed in accordance with the invention. 

[0027] FIG. 10 illustrates the change in acoustic reactance verses frequency along the 
x and y axis for the waveguide profile of FIG. 9. 

[0028] FIG. 11 illustrates the slope profile for each section of the waveguide 
illustrated in FIGs. 9 and 10 along the x and y-axis. 

[0029] FIG. 12 illustrates a depth verses height profile of another example waveguide 
designed in accordance with the invention. 

[0030] FIG. 13 illustrates the change in acoustic reactance verses frequency along the 
x and y axis for the waveguide profile of FIG. 12. 

[0031] FIG. 14 illustrates the slope profile for each section of the waveguide 
illustrated in FIGs. 12 and 13 along the x and y-axis. 

[0032] FIG. 15 illustrates a depth verses height profile of another example 
waveguide designed in accordance with the invention. 

[0033] FIG. 16 illustrates the change in acoustic reactance verses frequency along the 
x and y axis for the waveguide profile of FIG. 15. 

[0034] FIG. 17 illustrates the slope profile for each section of the waveguide 
illustrated in FIGs. 15 and 16 along the x and y-axis. 

[0035] FIG. 18 illustrates the acoustic frequency response of the waveguide shown in 
FIGs. 1 and 5 used with an electrical second order high pass filter, highlighting the 
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dispersion in the horizontal, vertical, and combination of horizontal and vertical 
directions. 

[0036] FIG. 19 illustrates the acoustic frequency response of the waveguide shown in 
FIGs. 1 and 5, highlighting the dispersion in the horizontal, vertical, and combination of 
horizontal and vertical directions. 

[0037] FIG. 20 is a flow diagram illustrating a design basis for a software program 
that performs according the methodology of the invention. 
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Detailed Description 



[0038] FIG. 1 illustrates a perspective view of a loudspeaker 100 utilizing an acoustic 
waveguide 102 designed according to the design method of the invention. As illustrated 
in FIG. 1, the loudspeaker system 100 has an acoustic waveguide 102 defined by a 
continuous three-dimensional surface. Defined at one end of the waveguide 102 is a 
throat 104 and at the opposing end, a mouth 106. Coupled to the throat 104 of the 
waveguide 100 is a transducer or driver 108. While FIG. 1, illustrates the loudspeaker 
driver 108 having a dome 110 diaphragm, loudspeakers using diaphragms of other shapes 
may also be used in connection with the invention. Further, the loudspeaker 100 in FIG. 
1, illustrates the waveguide 102 used in connection with a tweeter (generally 2kHz- 
20kHz); however, the waveguide 102 of the invention may be used in connection with 
specialized drivers for other dedicated parts of the audio frequency band, such as ultra- 
high frequency drivers (generally 10kHz-40kHz), midrange drivers (generally 200Hz- 
5kHz), and woofers (generally 20Hz-lkHz). 

[0039] FIG. 2 illustrates a cross-sectional view of the waveguide 102 and dome 
diaphragm 110 of the loudspeaker 100 of FIG. 1, taken along line A — A. As illustrated 
in FIG. 2, the throat 104 of the waveguide 102 is coupled to the diaphragm 110 of the 
driver. The waveguide 102 then flares outward from the throat 104 to the free end of the 
mouth 106 at an exponential flare rate m. 

[0040] FIG. 3 illustrates a flow diagram of an example implementation 300 of 
waveguide design methodology of the invention. As illustrated in FIG. 3, the initial step 
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302 of the invention involves establishing a set of performance metrics under which the 
waveguide is to be designed. In the described example implementation, the design metric 
under which the waveguide will be measured and designed is the minimum change in 
acoustic reactance. Although the design metric basis described in this example 
implementation is based upon the change in acoustic reactance, one skilled in the art will 
recognize that other design metrics, such as change in acoustic resistance, may be used in 
connection with the principles and theory of the invention to achieve substantially similar 
waveguide design profiles. 

[0041] Once the design metrics are established, an exponential waveguide profile 
with two or more different exponential slopes are then concatenated together 304. This is 
accomplished by first altering the slopes 306 of each section using the design metric. In 
this example implementation, the slopes are altered to sustain a constant change in 
acoustic reactance along the transitions section of the waveguide, from the throat to the 
mouth of the waveguide. Once the slopes 306 of each section are altered, the sections are 
then concatenated together using exponential functions based upon the desired depth and 
initial design radius of the given waveguide. Once the sections are concatenated 
together, the profile of concatenated exponential contours having modified slopes is then 
smoothed 308 based upon a polynomial order curve fit, producing a design key for ease 
in prototyping the waveguide. Steps 302-306 shall each be explained in further detail 
below. 
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[0042] FIG. 4 illustrates the transition in theoretical acoustic impedance along the 
wave front of a waveguide. As discussed in the background section, it is highly desirable 
to design a waveguide that can sustain a constant change in acoustic reactance or 
impedance along the transition of the waveguide from the throat to the mouth of the 
waveguide. Several know equations may be used to measure the change in acoustic 
reactance across a given waveguide and may be used as the design metric for the basis of 
the invention. 

[0043] To understand the equations defining acoustic impedance, it is first helpful to 
recognize several known theories associated with waveguides that may be considered to 
form the basis of the design metrics. The first equation of interest is: 

S = S T eT 

where St is the area at the throat, m is the flare rate along the length defined as x, and S is 
the area at the mouth of the waveguide. Further, steady state pressure is defined as: 

P{t) = P+e-"*e 21 } e- JQ)t 



where 



^ 2n tjj 27tf 



[0044] As for calculating or measuring the change in acoustic impedance across the 
waveguide from throat to mouth, it is known by those skilled in the art that acoustic 
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impedance is defined as unique components for low and high frequencies. For example, 
when the flare rate m is greater than 47i divided by the wavelength (m > 2k, low 
frequencies), the acoustic impedance is defined as: 



Rat =0 



Y - P ° C 

AT ~ ~7T~ 
O j 



( 



m 



2k \4k 2 \ 



and 



p 0 c = 406 mks ohms at 20° C and 1 0 5 newtons I m 2 ambient temperature 

[0045] Similarly, when the flare rate m equals 4n divided by the wavelength (m 
47i/ c / c where / c is the cutoff frequency) the acoustic impedance is defined as: 



R AT =Q 



y - p ° c 

* AT - ~ ~ 
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[0046] At this frequency, the acoustical impedance at all positions along the 
waveguide is reactive. As a result, no acoustical power will be transmitted below this 
frequency. 

[0047] Lastly, when the flare rate m is less than 4n divided by the wavelength, (m < 
2k, high frequencies), the acoustic impedance is defined as: 



[0048] As illustrated by FIG. 4, the acoustic impedance equation for high frequencies 
is generally applied near the throat 106 of the waveguide 102, where the waveguide 
interfaces with the diaphragm 110 of the driver. In contrast, the acoustic impedance 
equation for low frequencies is applied near the mouth 108 or open end of the waveguide 



[0049] The next step 304 of FIG. 3 is creating an exponential waveguide profile with 
two or more different exponential slopes concatenated together. To create this 




P 0 c 2 m 
2wS T 



102. 
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concatenated exponential waveguide profile, several input variable must first be 
provided, such as (i) the diameter or radius of the loudspeaker driver (or the initial radius 
or diameter of the throat of the waveguide); (ii) the initial slope of the waveguide along 
the x-axis (or major axis) and the initial slope of the waveguide along the y-axis (or 
minor axis); and (iii) the depth of the waveguide along the z-axis. FIG. 5 illustrates the 
x-axis, y-axis, z-axis and the radius r for an example waveguide for which the initial 
input variables may be obtained. 

[0050] Once the initial input variables are obtained, a waveguide having two or more 
concatenations may be created using the functions set forth below. As seen in FIG. 6, 
which illustrates a waveguide 102 divided into ten section, ml, m2, m3, m4, m5, m6, m7, 
m8, m9 and m 10, in the example embodiment, the waveguide 1 02 is divided into ten 
sections which are concatenated together as described in more detail below. The sections 
may be defined by sections of equal length along the depth of the waveguide. For 
example, if the depth is set at 1 inch, each section shall be 1/10 of an inch, or 0.10 inches. 
Although the example implementation divides the waveguide into ten sections, one 
skilled in the art will recognize that a waveguide profile may be obtained from two or 
more section that are concatenated together using the methodology of the invention as 
described in more detail below. 

[0051] The slope m of each section of the waveguide is derived by starting with the 
initial slope input along the x and y axis and modifying or updating the slope for each 
section along both axis such that the slope is optimal for a minimum change in acoustic 
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reactance a 1 1 he given f requency. T he o ptimal s lope for m inimum c hange i n a coustic 
reactance at a given frequency may be obtained from the derivate expressions of acoustic 
reactance at high and low frequencies. 

[0052] As previously discussed, for low frequencies, acoustic reactance is expressed 



as: 



X AT ~ 



Po c 



m l vi- 



lli \4k 7 



For high frequencies, acoustic reactance is expressed as: 



X AT ~ 



P 0 c m 
2wS T 



[0053] From these equations, the derivatives of the acoustic reactance with respect to 
low and high frequencies may be expressed, respectively, as the following functions: 



dx AT = Po c 

df S T 



1 



Srf J f mc 



f V 



Snf J 



-1 



4/sT 



, low frequency derivative 



dX AT _ p 0 c mc 



df 



S T 4nf 



, high frequency derivative 
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[0054] From these derivates, an optimization routine with respect to low frequencies 
for the slope m may be defined by the following: 



1 mc I mc 



-1=0 



Solving for m with respect to frequency we obtain an optimal slope for a minimum 
change in acoustic reactance at low frequencies as: 



m =^ = 4k 

Low frequencies 



[0055] While optimization for low frequencies is a helpful metric, it is advantageous 
to review an approach that considers the transition between low and high frequencies 
along the transition of the waveguide from the throat to the mouth of the waveguide. An 
alternate solution that captures a minimum change in the acoustic reactance between low 
and high frequencies along the throat of the waveguide may be defined as well. 
[0056] When equating the magnitude of the derivatives for respective low and high 
frequency expressions of acoustical reactance, we obtain: 



dX 



AT 



df 



PqC 



mc 

'mcY 



f 



-1 



mc 



p 0 c mc 
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[0057] Solving for m with respect to frequency, we obtain an optimal slope that 
minimizes discontinuities in acoustic reactance from low to high frequencies along the 
transition of the waveguide from throat to mouth: 



[0058] When considering the optimal slopes defined for a minimum in the derivative 
of the acoustic reactance, a specific slope update may be defined for each section of the 
waveguide. As such, the updates may be partitioned numerically into different regions. 
Optimal average updates can them be determined based upon the data g enerated from 
partitioning the waveguide in different regions. Smaller average slope updates may be 
used for the design of waveguides having shallower desirable depths. 
[0059] Tables 1,2 and 3 i llustrate t hree d ifferent p artitions o f w aveguides d ivided 
into ten sections having regions of particular interest between 1.5 kHz to 6 kHz. Each of 
the scaled slope updates in each of the different partitions of the tables are defined by the 
equations set forth below where the particular partitioning frequency divides the 
application of the equations: 



m 




16 2nf _ 16 



m 



= 4k and m 



16 2nf 16 




Low frequencies 



C 



Low to High frequencies 
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Slope Frequency(Hz) Slope Update 

Transition 



mlO(low) 1000 


0.136475364 


m9(low) 2000 


0.068237682 


m8(low) 3000 


0.045491788 


m7(low) 4000 


0.034118841 


Average Update 


0.071080919 


m6(low to high) 5000 


0.026428341 


m5(low to high) 6000 


0.022023617 


m4(low to high) 7000 


0.018877386 


m3(low to high) 8000 


0.016517713 


m2(low to high) 9000 


0.014682411 


ml(lowtohigh) 10000 


0.01321417 


Average Update 


0.01862394 


Table 1. Slope Update Partition (A). 


Slope Frequency(Hz) Slope Update 


Transition 




mlO(low) 1000 


0.136475364 


m9(low) 2000 


0.068237682 


m8(low) 3000 


0.045491788 


Average Update 


0.083401611 


m7(low to high) 4000 


0.033035426 


m6(low to high) 5000 


0.026428341 


m5(low to high) 6000 


0.022023617 


m4(low to high) 7000 


0.018877386 


m3(low to high) 8000 


0.016517713 


m2(low to high) 9000 


0.014682411 


m l (low to high) 10000 


0.01321417 


Average Update 


0.020682723 



Table 2. Slope Update Partition (B). 
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Slope 


Frequency(Hz) Slope Update 


Transition 






mlO(low) 


1000 


0 136475364 


m9(low) 


2000 


0 068237682 


m8(low) 


3000 


0.045491788 


m7(low) 


4000 


0.034118841 


m6(low) 


5000 


0.027295073 


m5(low) 


6000 


0.022745894 


Average Update 


0.05572744 


m4(low to high) 


7000 


0.018877386 


m3(low to high) 


8000 


0.016517713 


m2(low to high) 


9000 


0.014682411 


ml (low to high) 


10000 


0.01321417 


Average Update 


0.029025561 



Table 3. Slope Update Partition (C). 



[0060] As illustrated, the above Tables 1-3 are derived assuming a frequency range of 
1-10 kHz, commencing with the lowest frequency of 1 kHz applied to mlO or the mouth 
section. Thereafter, the frequencies are assigned in 1 kHz increments to each section, 
where the last section ml, or throat section, is assigned a frequency of 10 kHz. Table 1 
calculates the slope update for each of the ten sections (ml-mlO) by applying the 
equation for rate of flare at low frequency from a frequency range of 1 kHz to 4 kHz, 
then applies the equation for rate of flare from low to high frequency from 5 kHz to 10 
kHz. 

[0061] Table 2 calculates the slope update for each of the ten sections (ml-mlO) by 
applying the equation for rate of flare at low frequency from a frequency range of 1 kHz 
to 3 kHz then applies the equation for rate of flare from low to high frequency from 4 to 
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10 kHz. Similarly, Table 3 applies the equation for rate of flare at low frequency from a 
frequency range of 1 kHz to 6 kHz then applies the equation for rate of flare from low to 
high frequency from 7 to 10 kHz. While the above tables calculate the change in the rate 
of slope for particular regions of interest from 1.5 kHz to 6 kHz applying the design 
technique of using ten section to support contributions from 1-10 kHz, the above tables 
may be generated from other frequency regions of interest, such as 6 kHz to 12 kHz, 12 
to 20 kHz or 20 kHz to 40 kHz. 

[0062] The partitions outlined in Tables 1 , 2 and 3 may be used as a basis to update 
the slopes in update equations that may used in a waveguide design software code, two 
example commented code implementations of which may be found below. 
Slope Update Routine (A) 

% Composite exponential waveguide consisting of'i" sections, 
fori = 1:10, 
ifmlx<0.6 

mlx=mlx+0.0375*i; % Increased rate of flare. 
elseif mlx<0.7 

mlx=mlx+0.0575*i; % Considerably increased rate of flare. 
else 

mlx=mlx+0.0775*i; % Considerably increased rate of flare. 
end 

ifmly<0.6 

mly=mly+0.0275*i; % Increased rate of flare. 
elseif mly<0.8 

mly=mly+0.0675*i; % Considerably increased rate of flare. 
else 

mly=mly+0.0875*i; % Considerably increased rate of flare. 
end 
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Slope Update Routine (B) 

% Composite exponential waveguide consisting of n i n sections. 
fori = 1:10, 
ifmlx<0.6 

mlx=mlx+0.0175*i; % Increased rate of flare. 
elseifmlxO.7 

mlx=mlx+0.0275*i; % Considerably increased rate of flare. 
else 

mlx=mlx+0.0475*i; % Considerably increased rate of flare. 
end 

ifmly<0.6 

mly=mly+0.0275*i; % Increased rate of flare. 
elseif mly<0.8 

mly=mly+0.0475*i; % Considerably increased rate of flare. 
else 

mly=mly+0.0675*i; % Considerably increased rate of flare. 
end 

[0063] The slope update variables selected above may be based upon the average 
slope update obtained for each partition in the Tables above. The average update used is 
selected to be implementable, realizable and as close to an optimal solution as possible 
given the design parameters. Waveguides of large depths may accommodate greater 
rates of change since the transitions between the sections are larger. Thus, one skilled in 
the art may vary the update variables based upon the depth of the preferred waveguide 
design. Further, when the slope update routines are based upon elliptical waveguide 
designs, where the width of x-axis is typically longer than that of the y-axis, the rate of 
change in the slope may be more gradual along the x-axis than along the y-axis. When 
designing a circular waveguide, the rate of change may be equal along both axis, and thus 
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the routines for each axis may be identical, or as set forth above, may still vary producing 
a circular waveguide having elliptical dispersion patterns. 

[0064] Slope update routine A was developed to support alternate dispersion 
coverage for larger and deeper waveguides since the initial and cumulative rate of change 
is larger than Slope Update Routine B. Consequently, due to larger respective 
wavelengths, Slope Update Routine A would provide lower frequency response. Slope 
Update Routine B was used to implement a smaller and shallow waveguide with 
particular dispersion coverage. The two illustrations demonstrate the flexibility in the 
application of the methodology. Similarly, those skilled in the art could use a series of 
polynomial update functions to obtain a solution that provides performance in keeping 
with the design standard. 

[0065] For example, the above Slope Update Routine B was established for the 
purposes of creating a waveguide profile having a shallow depth of approximately 1 inch. 
Thus, the average slope updates of Table 3 were used as a basis for designing the Slope 
Update Routine B because of the rate of change in slope would produce a more optimal 
profile given the shallow depth of the waveguide. Thus, the average slope update rate for 
the low frequency of Table 3 was used as a guideline to define the two upper rates of 
change along the y-axis for the Slope Update Routine B (i.e. 0.05572744 falls between 
0.0475 and .0675). The lowest rate update along the y-axis was then used as the middle 
range update rate along the x-axis. The average slope update rate for the low to high 
frequencies of Table 3 was then used as a guideline to define the two upper rates of 
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change along the y-axis (i.e. 0.029025561 falls between 0.0275 and .0475). As will be 
demonstrated below, the above routine can also be used to design a waveguide having a 
mouth with a circular pattern, yet having elliptical acoustic dispersion characteristics. 
[0066] Using a slope update routine or formula implemented for an optimal solution 
given the design parameters, the slope of each section may be determined in a cumulative 
manner beginning with the initial slope input. For example, using the Slope Update 
Routine B, as set forth above, if the initial slope along the x-axis is 0.55 the initial slope 
ml will be updated to 0.56175 (0.55 + (0.0175*1)). For m2, the slope will be 0.59675 
(0.56175 + (0.0175*2)) and for m3 will be 0.64925 (0.59675 + (0.0175*3)). For m4, the 
slope will be 0.75925 (0.64925 + (0.0275*4)). Updates for sections m5-ml0 may follow 
cumulatively by the continued application of the Slope Update Routine B for the x-axis. 
The slopes for each section on the y-axis may be similarly calculated using the portion of 
the routine for cumulatively updating the slopes along the y-axis. 

[0067] Once all the slopes ml-mlO are established based upon the update formulas, 
the slopes may then be concatenated together based upon the initial radius of the desired 
waveguide profile at its throat and the depth of the desired waveguide profile design. 
One method for concatenating the sections together is to plot the radius (or height) of the 
waveguide along the x and y axis against the depth (z-axis) on a 100 x 10 matrix against 
the rate of flare m for each section defined by the updated slopes for each section. The 
first section ml in the matrix may be defined by 1:10, 1, section m2 by 1 1 :20, 2, section 
m3 is 21 :30, 3 and etc. Given a depth of 1 inch, the percentage of depth inches per each 



26 



U.S. Express Mail No. EU 124349650 US 
Filing Date: October 29, 2003 



PATENT 

Attorney Docket No. HI03025USU 
(P04009USU) 



point 1:100 on the matrix will be 0.01 inches. The height may then be determined along 
the matrix 1:100 (depth 1:100) based upon the established design metrics. For example, 
in this example implementation, the height at each 0.01 inches of depth from throat to 
mouth may be defined by the following equation: 

Outer_Radius(depth,subsection) = (Constant)*(Starting_radius*(exp(mlx*depth))) 
The concatenated sections ml-mlO can then be smoothed using a polynomial function 
order curve fit, using a order approximation that may set as an input variable, to create a 
waveguide profile, which may used as a design key for the waveguide profile. Both the 
concatenated sections and the smoothed concatenated section may be plotted on a matrix. 
Example of various plots of waveguide profiles are illustrated in FIGs 7-17, as explained 
in further detail below. 

[0068] FIG. 7 illustrates a depth vs. height profile of an example waveguide designed 
in accordance with the invention, using the Slope Update Routine A, set forth above. The 
initial parameters of the example profile of FIG. 7 include a throat radius of 0.55 inches, 
an initial slope on the x-axis of 0.55, and initial slope on the y-axis of 0.55, a depth of 1.0 
inch and a polynomial function order of 16. The concatenated updated slopes ml-mlO 
along the x-axis are illustrated by 702, while the concatenated updated slopes ml-mlO 
along the y-axis are illustrated by 704. The smoothed curve in accordance with the 
polynomial fit curve function is illustrated by 706 for the x-axis and by 708 for the y- 
axis. 
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[0069] FIG. 8 illustrates the change in acoustic reactance verses frequency along the 
x and y axis for the waveguide profile of FIG. 7. As illustrated by FIG. 8, the change in 
acoustic reactance along the x-axis 802 is the same and the change in acoustic reactance 
along the y-axis 804. The change in acoustic reactance along the x and y-axis may be 
calculated based upon the updated slopes for each section using the established design 
metrics for the change in acoustic reactance without the necessity of creating a prototype 
to test the performance of the waveguide. Although FIG. 7 illustrates a waveguide 
having an elliptical profile, the predicted changes in acoustic reactance along the x and y- 
axis of the elliptical waveguide are analogous to the performance of a waveguide having 
a circular design profile. 

[0070] FIG. 9 illustrates a depth verses height profile of an example waveguide 
designed in accordance with the invention, using Slope Update Routine B, as set forth 
above. The initial parameters of the example profile of FIG. 9 include a throat radius of 
0.50 inches, an initial slope on the x-axis of 0.50, and initial slope on the y-axis of 
0.1375, a depth of 1.0 inch and a polynomial function order of 16. The concatenated 
updated slopes ml-mlO along the x-axis are illustrated by 902, while the concatenated 
updated slopes ml-mlO along the y-axis are illustrated by 904. The smoothed curve in 
accordance with the polynomial fit curve function is illustrated by 906 for the x-axis and 
by 908 for the y-axis. 

[0071] FIG. 10 illustrates the change in acoustic reactance verses frequency along the 
x and y axis for the waveguide profile of FIG. 9. As illustrated by FIG. 10, the change in 
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acoustic reactance along the x-axis 1002 differs slightly from the change in acoustic 
reactance along the y-axis 1004. As before, the change in acoustic reactance along the x 
and y-axis may be calculated based upon the updated slopes for each section using the 
established design metrics for the change in acoustic reactance without the necessity of 
creating a prototype to test the performance of the waveguide. Although FIG. 9 
illustrates a waveguide having a circular profile, the predicted changes in acoustic 
reactance along the x and y-axis of the circular waveguide are analogous to the 
performance of a waveguide having an elliptical design profile 

[0072] FIG. 11 illustrates the slope profile for each section of the waveguide 
illustrated in FIGs. 9 and 10 along the x and y axis. The slope profile for the x-axis is 
represented by 1102 and illustrates the concatenated and smoothed slope profile of the 
updated slopes for sections ml-mlO along the x-axis. Similarly, the slope profile for the 
y-axis is represented by 1 104 and illustrates the concatenated and smoothed slope profile 
of the updated slopes for sections ml-mlO along the y-axis. 

[0073] FIG. 12 illustrates a depth verses height profile of an example waveguide 
designed in accordance with the invention, using Slope Update Profile B. The initial 
parameters of the example profile of FIG. 12 include a throat radius of 0.55 inches, an 
initial slope on the x-axis of 0.55, and initial slope on the y-axis of 0.1375, a depth of 1.0 
inch and a polynomial function order of 16. The concatenated updated slopes ml-mlO 
along the x-axis are illustrated by 1202, while the concatenated updated slopes ml-mlO 
along the y-axis are illustrated by 1204. The smoothed curve in accordance with the 
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polynomial fit curve function is illustrated by 1206 for the x-axis and by 1208 for the y- 
axis. 

[0074] FIG. 13 illustrates the change in acoustic reactance verses frequency along the 
x and y-axis for the waveguide profile of FIG. 12. As illustrated by FIG. 12, the change 
in acoustic reactance along the x-axis 1302 differs slightly from the change in acoustic 
reactance along the y-axis 1304. As before, the change in acoustic reactance along the x 
and y axis may be calculated based upon the updated slopies for each section using the 
established design metrics for the change in acoustic reactance without the necessity of 
creating a prototype to test the performance of the waveguide. Although FIG. 12 
illustrates a waveguide having a circular profile, the predicted changes in acoustic 
reactance along the x and y axis of the circular waveguide are analogous to the 
performance of a waveguide having an elliptical design profile. 

[0075] FIG. 14 illustrates the slope profile for each section of the waveguide 
illustrated in FIGs. 12 and 13 along the x and y axis. The slope profile for the x-axis is 
represented by 1402 and illustrates the concatenated and smoothed slope profile of the 
updated slopes for sections ml-mlO along the x-axis. Similarly, the slope profile for the 
y-axis is represented by 1404 and illustrates the concatenated and smoothed slope profile 
of the updated slopes for sections ml-mlO along the y-axis. 

[0076] FIG. 15 illustrates a depth verses height profile of an example waveguide 
designed in accordance with the invention, using Slope Update Routine A. The initial 
parameters of the example profile of FIG. 15 include a throat radius of 0.55 inches, an 
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initial slope on the x-axis of 0.55, and initial slope on the y-axis of 0.55, a depth of 1.0 
inch and a polynomial function order of 16. The concatenated updated slopes ml-mlO 
along the x-axis are illustrated by 1502, while the concatenated updated slopes ml-mlO 
along the y-axis are illustrated by 1504. The smoothed curve in accordance with the 
polynomial fit curve function is illustrated by 1506 for the x-axis and by 1508 for the y- 
axis. 

[0077] FIG. 16 illustrates the change in acoustic reactance verses frequency along the 
x and y axis for the waveguide profile of FIG. 15. As illustrated by FIG. 16, the change 
in acoustic reactance along the x-axis 1602 is the same and the change in acoustic 
reactance along the y-axis 1604. As before, the change in acoustic reactance along the x 
and y axis may be calculated based upon the updated slopes for each section using the 
established design metrics for the change in acoustic reactance. Although FIG. 16 
illustrates a waveguide having an elliptical profile, the predicted changes in acoustic 
reactance along the x and y-axis of the elliptical waveguide are analogous to the 
performance of a waveguide having a circular design profile. 

[0078] FIG. 17 illustrates the slope profile for each section of the waveguide 
illustrated in FIGs. 16 and 17 along the x and y axis. The slope profile for the x-axis is 
represented by 1702 and illustrates the concatenated and smoothed slope profile of the 
updated slopes for sections ml-mlO along the x-axis. Similarly, the slope profile for the 
y-axis is represented by 1704 and illustrates the concatenated and smoothed slope profile 
of the updated slopes for sections ml-mlO along the y-axis. 
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[0079] As illustrated above, by using the Slope Update Routines generated based 
upon the derived equations defining the optimal slope to minimizes discontinuities in 
acoustic reactance at low frequencies and from low to high frequencies, waveguides may 
be profiled to meet the performance standards of minimize the change in acoustic 
reactance along the waveguide at desired frequency. This is illustrated by measurements 
taken from a prototype fitting the design profile of the waveguide profile set forth in 
FIGs. 15-17. 

[0080] FIG. 18 illustrates the acoustic frequency response of the waveguide shown in 
FIGs. 1 and 5 used with an electrical second order high pass filter, highlighting the 
dispersion in the horizontal, vertical, and combination of horizontal and vertical 
directions. The line identified as 1802 represents the on-axis response, which is generally 
defined as the direct radiating contribution. Line 1804 represents the listening window, 
which is representative of dispersion in nominal listening conditions. Line 1806 
represents first reflection, which is the dispersion with annular interaction with listening 
room walls. Line 1808 represents sound power, which is dispersion with contribution of 
the waveguide in 360 degrees. Line 1810 is the directivity of sound power, while line 
1812 is the directivity of first reflection. The directivity of sound power is the on-axis 
response s ubtracted f rom t he s ound p ower d ispersion a nd d efines u niformity o f s ound 
power dispersion or overall dispersion referenced to on-axis contribution. Directivity of 
first reflection is the on-axis response subtracted from the first reflection and defines the 
uniformity of the first reflection dispersion reference to on-axis contribution. 
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[0081] FIG. 19 illustrates the acoustic frequency response of the waveguide shown in 
FIGs. 1 and 5, highlighting the dispersion in the horizontal, vertical, and combination of 
horizontal and vertical directions. In this FIG. 19, the on-axis response 1802 of FIG. 18, 
listening window 1804, first reflection, 1806, and sound power 1808 are overlaid on top 
of one another, represented by 1902, to demonstrate the similarities of the dispersions 
patterns of the waveguide under the various responses, in the frequency range of interest, 
which in this illustration is 1kHz to 10kHz. Similarly, lines 1810 and 1812, the 
directivity of sound power and directivity of first reflection, respectively, are overlaid, as 
illustrated by 1904, to demonstrate the similarities in dispersion patterns for these 
measurement in the frequency range of interest. As illustrated by FIGs. 18 and 19, the 
waveguide designed in accordance with the invention provides uniform coverage over a 
wide range of frequencies. 

[0082] Further, as illustrate above, the Slope Update Routines A and B may be used 
to design waveguides having elliptical cross-sectional areas that produce circular 
dispersion patterns (i.e. an elliptical waveguide that produces the same horizontal and 
vertical dispersion patterns from low to high frequencies, which, by way of example, may 
be considered between 1kHz to 10kHz). Conversely, the design methodology allows for 
the profiling of waveguides having circular cross-sectional areas but that provide 
elliptical dispersion patterns (i.e. a circular waveguide that produces different horizontal 
and vertical dispersion patterns from low to high frequencies, which, maybe considered 
in range between 1 kHz to 1 0kHz). T o design an elliptical w aveguide h aving circular 
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dispersion characteristic, the initials slopes along the x and y-axis are substantially the 
same. In contrast, to design a circular waveguide having elliptical dispersion 
characteristics, the slopes along the x and y-axis may initially differ. 
[0083] While the above described profiles can be generated by hand calculations and 
plotted from raw data, the above described methodology can be embodied in a software 
program that will calculate and the plot output data and smooth the curves based upon a 
polynomial order function. The software can also create output files containing the raw 
data as well as the plotted curves. The process may be performed by hardware or 
software and may be designed within known software programs, such as Matlab®, a 
software program sold under the registered trademark by The MathWorks, Inc., or may 
be designed as a standalone executable software program. 

[0084] FIG. 20 illustrates a flow diagram that may be used as a design basis for a 
basic software program that performs according the methodology of the invention. One 
skilled in the art will recognize that the ordering of functions set forth in FIG. 20 may 
vary. For example, the initial input variable may be given before the design metric is 
chosen. F urther, t he p rogram m ay be d esigned t o give t he u ser t he c hoice o f s everal 
design m etrics. A dditionally, t he c ertain v ariables m ay b e p redetermined, s uch a s t he 
number of sections of the waveguide and the order number for the polynomial fit curve, 
or the variable may be input by the user. 

[0085] As illustrated by the FIG. 20, the program may first ask for design input 
variables, such as initial t hroat r adius, initial slope along the x and y-axis, waveguide 



34 



U.S. Express Mail No. EU 124349650 US 
Filing Date: October 29, 2003 



PATENT 

Attorney Docket No. HI03025USU 
(P04009USU) 



depth, number of sections of the waveguide, the design metric and/or the polynomial 
order fit number 2002. Many of these design variables may, however, be predetermined, 
such as the depth of the waveguide, the number of sections of the waveguide, the 
polynomial order curve fit number and the design metric. 

[0086] Once the input variable are collected, the program may then update the 
sections of the waveguide based upon the design metric 2004. If the design metric is 
predetermined, such as the change in acoustic reactance, the program may be design with 
predeveloped software routines, such that the Slope Update Routines A and B set forth 
above. However, these routines may be developed by the program based upon the input 
variables and then applied to update the slopes of the sections of the waveguide. 
Additionally, the software could be designed using other software routines for different 
design metrics, or may give the user the option of selecting different design metrics for 
the design of the waveguide, such as change in acoustic resistance, change in acoustic 
resistance, minimum change in acoustic resistance, minimum change in acoustic 
reactance, and then applying different slope updates for the selected design metric. 
[0087] Once the slopes for each section have be defined or updated by the application 
of the design metric, the sections can then be concatenated together using known 
equations f or d etermining t he p rofile o f a w aveguide given t he s lopes of e ach s ection 
2006. The concatenated sections can then be smoothed using a polynomial order fit 
curve or other similar method 2008 to create a design key for the waveguide. 
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[0088] The profile can then be validated by calculating the performance of the 
waveguide based upon the design metric 2010. Although not shown on FIG. 20, design 
iterations may be made to the profile if the performance calculations are undesirable or 
need improvement. The waveguide profiles and performance standards can be plotted 
and output files, such as excel spreadsheets, data files, text files or tables, can be 
generated from the raw and smoothed data 2012. One skilled in the art will recognize 
that t he d esign i terations m ay b e m ade t o a s oftware p rogram o f t his t ype t o add a nd 
remove features, to make the program more user friendly, to provide user options or to 
use the program to calculate or determine the slope update routines depending upon the 
design metric. 

[0089] The process described above may be performed by hardware or software and 
may be designed within known software programs, such as Matlab, or may be designed 
as a standalone executable software program. If the process is performed by software, 
the software may reside in software memory (not shown) in the controller 1012, memory 
device 1014, Call Processor 1006, GPS module 308, or a removable memory medium. 
The s oftware i n m emory may i nclude an o rdered 1 isting o f e xecutable instructions for 
implementing logical functions (i.e., "logic" that may be implement either in digital form 
such as digital circuitry or source code or in analog form such as analog circuitry or an 
analog source such an analog electrical, sound or video signal), may selectively be 
embodied in any computer-readable (or signal-bearing) medium for use by or in 
connection with an instruction execution system, apparatus, or device, such as a 
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computer-based system, processor-containing system, or other system that may 
selectively fetch the instructions from the instruction execution system, apparatus, or 
device and execute the instructions. In the context of this document, a "computer- 
readable medium" and/or "signal-bearing medium" is any means that may contain, store, 
communicate, propagate, or transport the program for use by or in connection with the 
instruction execution system, apparatus, or device. A signal-bearing medium 
encompasses a computer-readable medium. The computer readable medium may 
selectively be, for example but not limited to, an electronic, magnetic, optical, 
electromagnetic, infrared, or semiconductor system, apparatus, device, or propagation 
medium. More specific examples "a non-exhaustive list" of the computer-readable 
medium would include the following: an electrical connection "electronic" having one or 
more wires, a portable computer diskette (magnetic), a RAM (electronic), a read-only 
memory "ROM" (electronic), an erasable programmable read-only memory (EPROM or 
Flash memory) (electronic), an optical fiber (optical), and a portable compact disc read- 
only memory "CDROM" (optical). Note that the computer-readable medium may even 
be paper or another suitable medium upon which the program is printed, as the program 
can be electronically captured, via for instance optical scanning of the paper or other 
medium, then compiled, interpreted or otherwise processed in a suitable manner if 
necessary, and then stored in a computer memory. 

[0090] The software may be processed by a processor such as general purpose 
microprocessor, application specific processor ("ASP"), digital signal processor ("DSP"), 
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application specific integrated circuit ("ASIC"), and/or reduced instructions set integrated 
circuit ("RISC") processor. 

[0091] While various embodiments of the invention have been described, it will be 
apparent to those of ordinary skill in the art that many more embodiments and 
implementations are possible within the scope of this invention. For example, the same 
principles used to design waveguides, as described herein, may be found used in radar 
and communication applications using analogous partitions, concatenations, and design 
metrics. Further, the waveguides may be the diaphragms of the loudspeakers. The design 
approach may also be applied to the design of port tube profiles found in loudspeaker 
systems, as well as, waveguides. The design approach in connection with the port tubes 
would require the application of appropriate functions for the desired port tube flare rates 
that are concatenated contributions of respective metrics that overall, increase the useable 
headroom of the port tubes in loudspeakers. Accordingly, the invention is not to be 
restricted except in light of the attached claims and their equivalents. 
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